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ABSTRACT 



A global, six-level, primitive eguatian model was numer- 
ically integrated to examine the interaction of a mid-lati- 
tude baroclinic wave with several different mountain-valley 
configurations. The purpose was to determine the contribu- 
tion of a well-developed baroclinic wave to the initiation 
of a monsoon surge. It was found that these waves can ini- 
tiate a surge. Such surges are wealc and limited in extent 
but they do exhibit the characteristics of observed monsoon 
surges. The latitudinal location of the east-west mountain 
range has a large impact on the strength of the disturbance 
and thus on the initiation of the surge. However, the size 
of the valley is not an important consideration. The surge 
is strongest at the lowest levels and dissipates rapidly 
with height. Overall results indicate surge forcing by the 
baroclinic wave is possible but other forcing mechanisms are 
necessary in crder to simulate the stronger, more extensive 
surges observed during the winter monsoon. 
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I. 



INTROD UCTION 



In recent years meteorologists have given considerable 
attention to the Asian winter monsoon. Described by Chang 
and Lau (1980) as "one of the most energetic heat engines in 
driving the earth's atmosphere", the winter monsoon is the 
most distinctive component of the Northern Hemispaere winter 
circulation. Its fundamental, large scale structure has 
long been observed, however, much about the monsoon remains 
unknown. The role of the monsoon in the interaction of the 
mid-latitude and tropical regions, the variability of the 
monsoon and the forcing cf small scale monsoonal variations 
are all areas of active research. 

A monsoon may be defined as a three dimensional, plane- 
tary scale wind regime that reverses direction seasonally. 
In the case of the monsoon, the establishment of the summer 
cyclone and winter anticyclone is governed by land-sea dif- 
ferences and the varying position of the sun. Essentially, 
it is the seasonal replacement of one persistant circulation 
system by an equally persistant but opposite system. The 
development of these circulations closely links the monsoon, 
particularily that of east Asia, to regional topography. 
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The long, clear nights of the Siberian winter allow 
strong radiaticnal cooling and the subsequent formation of a 
shallow, cold layer of surface air. Persistent, cold air 
advection throughout the troposphere further contributes to 
the formation of a surface high. The winter monsoon, which 
dominates the entire Asian continent, is characterized by 
this intense, cold polar anticyclone. The resulting north- 
northeasterly surface flow extends over all of eastern Asia 
and southward across the equator. Tha Himalayan Mountains 
restrict the southward movement of tha cold lower tropos- 
pheric air and provide an effective barrier to the mixing of 
the cold continental airmass with the warm tropical airmass 
of southern Asia. As the anticyclone becomes established 
over Siberia, a surface cyclone develops over the "maritime 
continent" region of Indonesia and Malaysia as a major plan- 
etary scale convective area is shifted eastward from its 
summer position near India (Ramage, 1971) . 

An analysis of the monthly mean velocity potential field 
at 200 mb for December 1974 by Chang and Lau (1980) shows a 
broad region of divergence over the maritime continent. 
This is indicative of deep rising motions and large scale 
overturning of the atmosphere. Their analysis also shows a 
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region of maximum convergence over nortiscn China in the 
vicinity of rhe Siberian anticyclone. (Sote char for large 
scale flow the divergence pattern is appr oximarel y given by 
the velocity potential.) Upper-level convergence is also 
evident in the equatorial west Pacific 3cean and over the 
east African coast. Although much weaker, these latter two 
areas indicate the existence of a '^alJcer circulation. 
Ascending over the maritime continent, one branch flows 
westward across the Indian Ocean and Arabian Sea, while 
another branch flows eastward to descend over the western 
Pacific Ocean. 

Thus the winter tropospheric circulation in this region 
is dominated by a strong local Hadley cell; a thermally 
direct circulation involving descent in the cold air over 
the northern anticyclone, southward surface flow, ascent of 
warm air over the maritime continent and northward return 
flow at 200 mb. 

The large-scale components of the monsoon are consistent 
and identifiable from year to year, but they are not static. 
The primary features not only vary in strength and location 
on an annual time scale but there are small time scale vari- 
ations in the velocity, temperature and pressure fields as 
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w9ll. One particular event that can be linicad to short rime 
scale variations is the mcnsoon surge. Also known as "cold 
surges" rhese events occur a- intervals of several days to 

two weeks during the winter monsoon (Lim and Chang, 1981). 

The characteristics and effects of the monsoon surge 
have bean discussed in several observational studies at both 
the synoptic and planetary scales by Danielson and Ho 

(1 969), Chang et ^ ( 1979), Murakami (1 979), Chang and Lau 

(1980,1982) and Chang et al (1983). The following sequence 
of events for a monsoon surge is taken from these studies. 

Prior to the occurrence of a monsoon surge there is 
intense, deep cold air advection over northern Asia accompa- 
nied by the passage of a mid-tropospheric trough that inten- 
sifies rapidly over northern Japan. The resulting increase 
in baroclinicity is associated with anticyclogenesis over 
Siberia and cyclogenesis in the strong baroclinic zone off 
the east China coast. The east Asian Hadley circulation is 
strengthened through increased subsidence as the anticyclone 
intensifies. The east Asian subtropical jet accelerates 
near Japan as the intensification of the Hadley circulation 
increases upper-level ageostrophic flow. The west Asian jet 
streak varies inversely with the east Asian jet and so 
decreases in strength over Afghanistan and Pakistan. 
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The monsoon surge is initiated as the upper-level trough- 
moves eastward and the pressure gradient across the east 
China coast rightens. The normal geostrophic balance of rhe 
anticyclonic flow is broken since the pressure gradient 
force is now greater than the Coriolis force and a rapid 
cross-iso baric flow from high to low pressure begins. This 
accelerating flow towards lower pressures at lower tropos- 
pheric levels is the monsoon surge. Following the surge, an 
increase in the northeasterly flow can be traced to the 
equatorial regions within 12-24 h. Within one to two days 

convection in the equatorial region intensifies and sustains 
the already enhanced Hadley circulation. The upper-level 
outflow from the equatorial South China Sea accelerates not 
only northward, but within one day of the surge accelerates 
east and west along the equator, thus strengthening the 
Walker circulation. The surge ends as the mid-latitude 
trough/ridge pattern moves sufficiently eastward so that 
subsidence no longer occurs over the anticyclone. 

Lim and Chang (1981) and Chang, Chan and Hillard (1983) 
suggest that the monsoon surge is characterized by two 
stages separated by an interval of a few hours to one day, 
depending on the location of the observing station. The 
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first stage is the pressure surge, or leading edge of the 
downward accelerating air, with a propagation speed charac- 
teristic of gravity waves. The second stage, which aoves at 
advective time scales, is due to a frontal passage and is 
characterized by a sharp decrease in the surface dewpoint. 

Few studies have been dona to determine the cause or 
nature of the forcing that precipitates a monsoon surge. 
Lim and Chang (1981) used the linearized shallow-water equa- 
tions on an equatorial beta-plane in an analytical study of 
monsoon response to a transient mass source that simulated 
intense anticyclogenesis. Their results showed that such 
mid- latitude forcing excites mainly Rossby waves after an 
initial period of gravity wave type motions associated with 
strong northerly winds. They suggested that pressure forc- 
ing, using equatorial beta-'plane dynamics, can explain the 
large features of synoptic scale motions in the eastern mon- 
soon region. However, they do not consider other effects 
such as orography or a variable mean zonal wind. Nor does 
their method account for the shallow vertical structure of 
monsoon surges or upper tropospheric response to the surge. 

The purpose of this study is to examine the interaction 
of topography, specifically a high east-west mountain range. 
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with a mid-iatituds baroclinic wave. Th= purpose will be zo 
determine the contribution of the baroclinic wave ro the 
initiation of monsoon surges. A global, six-level primitive 
equation model will be used to simulate the meteorological 
conditions prior and .subsequent to surge initiation. Ini- 
tial conditions will include a strong baroclinic jet north 
of the mountain range and a small disturbance in the mean 
current; the disturbance will grow through baroclinic insta- 
bility. Several mo un tain -valley configurations will be used 
to study the surge and numerical solutions will be compared 
with observed monsoonal surges. 
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II 



. DSS CRI PTION OF MODEL 

The numerical model used in these experiments is a ver- 
sion of the UCLA potential enstrophy conserving, general 
circulation model for an adiabatic, dry and frictionless 
atmosphere. The governing equations for this model are the 
primitive equations; the horizontal momentum equation, 
thermodynamic energy equation, surface pressure tendency 
(continuity) equation, and the hydrostatic equation. These 
are given in Appendix A. The design and dynamic processes 
of the model are fully discussed by Ara!cawa and Lamb (1977, 
1981) . A summary of the version used for this experiment is 
presented here. 

A. THE HOEI2CNTAL GRID 

The prognostic variables, wind(U,V) , temperarure (T) , 
and terrain pressure (n), are staggered horizontally accord- 
ing to Arakawa's Scheme 'C (Fig. 1). This grid arrange- 
ment, although using more computer time than an unstaggered 
grid, gives a better structure for the shorrer waves that 
arise due to geostrophic adjustment, nonlinear interaction 
or mountain effects. Resolution is 2.5° long by 4.0® lar 
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using a 36 x 46 spharical gril. Tha conpu tational domain is 
a 90*^ sector of the globe wi-h cyclic continuity imposed at 
the east and west boundaries. 

In spherical coordinates the poles are singular points. 
Since velocity components cannot be defined at the poles, 
they are taken as n-points. Due to the convergence of the 
meridians towards the poles a very short time step would 
normally be required to insure computational stability. To 
allow a longer time step, longitudinal averaging is done for 
selected terms in the equations that govern gravity wave 
propagation. To do this the zonal pressure gradient and 
zonal mass flux are expanded into a Fourier series around 
each latitude circle, then the amplitude of each wave compo- 
nent is reduced by the factor S, This technique smooths the 
gradients of the meteorological fields rather than their 
Fourier expansions. The smoothing operator S is applied 
only at higher latitudes (where S<1) and is given by; 

_ _ ^ cos<t> 

A(|i sin(nAA/2) 

(Symbol definitions are given in Appendix 3> 
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B. THE VERTICAL GRID 

The equations used in this model are written using the 
vertical sigma coordinate: 

p - p 

a = t where II = P “ P^ 

"T— = 

Ps is the surface pressure and is the pressure at the 

tr opopause . 

The vertical domain extends from the earth's surface to 
200 mb. Variables are staggered vertically (Fig. 2) in six 
layers equally spaced in sigma. The lower boundary, cr = 1, 
follows topography, where the surface pressure is defined. 
At both the upper and lower boundaries vertical velocity, a, 
vanishes. The prognostic variables (u,v,I,) are computed at 
mid-layer, while gaopoten tial , g, vertical velocity, a, and 
terrain pressure, n are computed diagnostically on sigma 
le vels. 

C. FINITE DIFFERENCING SCHEMES 

The time differencing methods used in this model are the 
fourth-order space differencing scheme and the Matsuno 
scheme (also called the Euler-backward scheme) . One Matsuno 
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time step is performed followed by four centered time steps. 
The Matsuno scheme is used to eliminate the computational 
mode of the centered scheme and because it selectively damps 
high frequency waves. The time step used here is 360 s. 
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Fi gure 1 . 




Horizontal Distribution of Variables on a 
Spherical Grid. 
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III. 



y.Q DEL INITI ALIZ ATION 

Th5 model uses an analytical initialization scheme 
developed by Major J. L. Ha yes, ^ OSAF. A small perturba- 
tion is superimposed on a horizontally and vertically vary- 
ing mean current, then allowed to grow through baroclinic 
instability. The scheme was derived specifically for flow 
over a flat earth; initialization over high terrain requires 
an adjustment to the scheme and is discussed separately. 

A. INITIALIZATION OF THE MEAN FLOW 

A baroclinic Bicklay jer is situated at 45°N, where it 
reaches a maximum velocity of 46 m/s at 200 mb. The initial 
mean surface flow has the same latitudinal variation as the 
jet; it reaches a maximum of 5 m/s at 45‘^N and decreases to 
zero at the poles. The v-component of the wind is initially 
set to zero. The surface wind is described by the Bicklay 
jet profile: 

U (<b) = U sech^ Cy ) 

S lOL 

V (4>) =0 

s 



‘ Personal communication, present affilitation Detachment 1, 
HQ AWS. 
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The upper-level wiad has the same latitudinal profile: 

The last term cn the right produces the vertical shear. The 
initial wind at any latitude and level is calculated by 
adding the upper-level and surface winds. 

Since the mean flow is in geosrrophic balance, the sur- 
face pressure is calculated by integrating the geostrophic 
wind equation: 

?(<}») = P (ijio) 

s s'' 



exp I- 



(■ 



tanh(y ^ J 

RTqY 






= sech^ (y 






’ln(P^/P) 

ln(P /P ) 
s m ' 



= 0 



The value of T, , 288K, is taken from the NACA standard 

atmosphere using a standard surface pressure of 1013.25 mb. 

Temperature, except at the surface, is determined by 
integrating the geostrophic thermal wind equation: 



T(p,(J)) 



T(p,<{)o) 



0 a 2Q 
m 

R In CP/P ) 
m 



I 



sinij) sech^ (<{)-<{)o) d({) 

<i>0 



Integrals are computed using the Simpson Rule approximation. 
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The vertically varying term T(p,?>.) is obtained from tae 
NACA standard atmosphere: 



T(p,<^o ) 



= T(po,<J>o) ~ 



where ; 




The constants are defined in Appendix B. 



3. INITIALIZATION 07 3H HIGH TERRAIN 

The mountains in this model are intended to be a very 
simple representation of the Himalayan Range. The mountains 
extend east to west across the grid section with a maximum 
elevation of 3000 m. The longitudinal span is 16® or 
approximately 1600 km. 

The mountains were initialized at full height without 
any growth period. The equations used to produce the moun->- 
tains are relatively simple and produce a knife-edged range: 



$ = = 



.25 at = 2A(j) 

. 75 at ] I = 2A(J) 

at .}) = 



The immediate initialization of the terrain does lengthen 
the geostrophic adjustment time. However, the model is run 
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for 72 h prior to creating the valley and it is reasonable 
to assume the longer adjustment period will have no serious 
effect on the attempted simulation of the monsoon surge. 
Using the mean flow initialization scheme without a mountain 
growth period resulted in a large scale, north-south oscil- 
lation in the pressure field. There was little evidence of 
damping during a 96-h model run. To solve this problem the 
surface pressures at mountain grid points were adjusted for 
height above mean sea-level. Using the hydrostatic equation 
and the adiabatic lapse rate for the UACA standard atmos- 
phere, the following equation was developed for a new (moun- 
tain) surface pressure in terms of the old (flat earth) 
surface pressure: 



The new surface pressure is used to re-calculate wind and 
temperature values at each mountain grid point. No oscilla- 
tion was evident when the adjusted values were used in the 
initialization. 

As the monsoon surge is expected to occur almost immedi- 
ately following the creation of the valley, the geostrophic 



P (new) 
s 




) 
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3^-fl 



adjustment nime cannot be neglected an this point. There- 
fore, the valley is “built" over a 12-h period no minimize 
the time required for geostrophic adjusnmenn. The menhod 
used is based on the scheme developed by Hayes (previously 
cited) for raising terrain height during a model run. To 
build the valley however, the method is reversed: 

* = Cl fl - C2 (T)G(X)'| 



The function C,(t) is defined: 



C2(t) = 



• 2A \ 

sxn^j— ) 

V= ■f-\J 



The funcnion G(X) is defined: 



G(X) = 



X-Xi-2AX 

-2AX 

X 7 + 2 AX-X 

2AX 

0 



i f 

and t>t. 
1 f 

Xi<X<X2 

'(X-2AX)^X^Xi 

X2^X^(X2+2AX) 

elsewhere 



X, and X^ are the west and east edges of the 
tively. Where no valley is built, G (X) - 

walls are sloped over two grid-lengths (2^X) 
floor is at mean sea- level (0 m) . 



valley respec- 
1 . The valley 
and nhe valley 
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c. 



TH2 PERTURBATION 



The perturbation is a standard barorropic 
It consists of a weak longitudinal sine wavs 
latitudinal amplitude at 45°N. The vertical 
one-tenth the earth's radius. The perturbation 



disturbance . 
with maximum 
amplitude is 
is described 



by: 



$* = fo (a sin(nX) sin^(2(Ji)) 



The following perturbation quantities are derived using the 
hydrostatic equation and geostrophic wind equation: 



p’ 

u* 

V’ 

T’ 



" 'ot " - H’(a sin(nl) Sin2(2(|))) 



= ■ 2sin(44>)) 



a£« cos 6 ij) 
0 



6^ * _ 1 i^nA cos(nX) sin^(2<j>)^ 



a coS({) 



The total fields 
quantities to the 
be integrated in 
values. 



are calculated by adding the perturbation 
aean current. The total fields will then 
time to arrive at the required forecast 
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IV. 



ANALYSIS OF RES NLIS 

A. EXPERIMENTS 1 S 2; CONTROL RUNS 

Twc ccntrcl runs were made for comparison wirh subse- 
quent model runs. The first run uses a flat earth while the 
second includes a 3000-m east-west mountain range centered 
at 30®N. Both runs are initialized with the barcclinically 
unstable mean flow and the barotropic disturbance discussed 
in Chapter 3. Model times given are based on the time (h) 
at which the disturbance is introduced into the mean flow. 

The first run follows the development of the disturbance 
independent of terrain affects. Sea-level pressure fields 
for 60, 12 , 84 and 96 h are shown schemetically in Figs. 3 

and 4. Potential temperature fields for the same times are 
shown in Figs. 5 and 6. Since a pie-shaped section of the 
globe is represented in cartesian coordinates there is some 
distortion of the fields. This is primarily evident in the 
tight north-south pressure gradients waioh, for comparable 
wind velocities, appear weaker in the east-west direction. 

The disturbance is strictly limited to the mid-latitudes 
with a very flat pressure gradient in both the polar and the 
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tropical regions. The disturbance moves eastward at a con- 
stant velocity of approximately 8.5° every twelve hours. 
Wave development is realistic wirh a 995 mb low and a 1023 
mb high developing after 78 h. 

Development of the disturbance in the wind and tempera- 
ture fields closely parallels that of the pressure field, 
although the temperature and potential temperature data dis- 
play far more short wave noise. As no surface frictional 
dissipation, diffusion, or smoothing techniques are used in 
the model, some noise in the results is expected. The large 
scale temperature features remain discernable throughout the 
six day (144 h) run but the the small scale features are 
masked by noise after 126 h, especially at higher levels. 
(Temperature and wind velocities are plotted on sigma sur- 
faces; 'higher levels' refers to increased height, or lower 
sigma values. ) 

After 96 h par ti cularily strong cyclogenesis begins to 
occur. At 114 h the cyclones intensify rapidly and the 
ridge is drastically weakened. While diffusion or smoothing 
applied to the fields would decrease the strong gradients 
caused by cyclogenesis, they would also weaken the gradients 
of pressure and temperature which could contribute to possi- 
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bl9 monsoon surges. For this reason neither method is used. 
A linear drag term was applied to the winds in the lowesr 
level, however it had no discernable effect on rhe strong 
cyclogenesis. 

The mountain range used in second control run is the 
basis for all subsequent terrain configurations. Fig. 7 
gives a three-dimensional representarion of rhe range. 
(Mote the vertical scale is greatly exaggerated.) The data 
from this run are used to determine the rime and location at 
which the valley is to be created. The requirements are: 
(1) a high pressure center to the northeast (to be consis- 
tent with the Siberian anticyclone) and (2) a sea-level 
pressure difference of 25 mb between the high and low pres- 
sure centers. The latter requirement ensures an adequate 
northwesterly currant to initiate a monsoonal surge. 

Sea-level pressure and potential temperature fields for 
60, 72, 84 and 96 h are shown in Figs. 8 through 11. East- 
ward movement of the disturbance is again 8.5® per twelve 
hours. The most noticeable difference between this experi- 
ment and the first is a series of weak highs and lows that 
have developed along the northern slope of the mountains. 
These disturbances migrate eastward at a speed comparable to 
the mid-latitude wave. 
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Devalopnient of rhs wave is slower thm in the previous 
case. At 78 h a 999 mb low and a 1019 mb high have devel- 
oped, giving a pressure difference of 20 mb compared to 28 

mb earlier. This represents an overall weakening of the 
disturbance and is not due to the weakening of either the 
high or low centers individually. It appears the mountain 
range may have a stabilizing effect on the growing distur- 
bance. The required 25 mb pressure difference is reached at 
84 h. Since a twelve hour period is required to create the 
valley, growth will be initiated at 72 h model time in sub- 
sequent runs. Since the valley must be located east of the 

anticyclone in order to taka advantage of the associated 

northwest flow, the western edge of tne valley is placed 
directly south of the leading portion or northwesterly winds 
at 84 h. 

Development of the disturbance in the potential temperar 
ture and wind fields is very similar to the first run. 
There is a maximum in potential temperature over the moun- 
tains but this is simply due to the increased elevation. 

The small scale noise apparent in the potential tempera- 
ture fields of the first run also appear in this experiment 
although at a later time, possibly due to the much slower 
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development of the cyclone-anticyclone system. The s~rong 
cyclogenesis and intense pressure gradients also appear 
lexer in the second model run. Wave development is approxi- 
mately six hours slower when the mountains are included and 
this corresponds to the time lag observed in the cyclogene- 
sis and noise growth. 

It should be noted that potential temperature is used in 
the analysis of model results since its fields will not 
include effects of the adiabatic processes occurring in flow 
over mountains. 

B. EXPERIMENT 3 

The initialization for this experiment, and all subse- 
quent experiments, is provided by the same analytic scheme 
used in the control runs. Terrain consists of an east-west 
mountain range centered at 30<^N. The valley, built between 
72 and 84 h, is twelve gridlengths wide or nearly one-third 
the length of the mountain range (Fig. 12) . The critical 
time period for the occurrence of a surge will be the first 
24 h after the valley is completed (84-108 h) . It is during 
this period that the anticyclone is located northwest of the 
valley. 
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Development of the disturbance through 72 h is described 
by control run f2. The sea-level pressure field a.z 84 h 
shows a trough extending across the valley and along the 
northern slope of the eastern range. At the same time a 
small, high pressure ridge can be seen in the northwest 
region of the valley. For the next twelve hours -he ridge 
continues to build southward until, at 96 h, it extends the 
length of the western valley wall splitting the trough into 
two separate low centers. The pressure gradient in the 
tropical region south of the valley has also increased (Fig. 
13) . As the ridge continues to expand southward, the east- 
ern portion of the trough develops a very strong, tight 
pressure gradient in the valley and along the northern moun- 
tain slopes (Fig. 14). The ridge continues to strengthen 
and remains located along the western wall through 114 h. 
After this time the anticyclone moves eastward to a position 
due north of the valley and the ridge begins to weaken. 

Also noted in the results is an oscillation in the 
intensity of the pressure field. It begins when valley 
growth is initiated and has a period of approximately 24 h. 
The observed oscillation is possibly the result of gravity 
waves set up by the relatively sudden (12 h) creation of the 
valley. 
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The second stage of the monsoon surge was earlier 
described as being characterized by a sharp drop in surface 
dewpoint. Since this is a dry model, and as observational 
studies (Chang et al, 1979, Chu, 1978) show rhis stage can 
also be associated with a drop in surface temperatures, the 
potential temperature fields are analyzed for possible surge 
characreristics. (Beyond this point ' teaperarure • refers to 
potential temperature unless otherwise stared.) 

In agreement with the development of the surface pres- 
sure ridge, a weak tongue of cold air has formed along the 
western valley wall at 84 h. However, no similar gradient 

exists along the eastern wall (Fig. 15). 

Twelve hours after the valley is completed, the cold air 
has reached the southern valley boundary. Minimum tempera- 
ture at 24®N has dropped from 310K at 34 h to 303K at 96 h. 
A tight temperature gradient exists along the western side 
of the cold air intrusion and also along the northern moun- 
tain slopes where temperatures have remained relatively con- 
stant. The cold air tongue can easily by seen at the two 
lowest levels (Fig. 16) . At the mid-levels a very slight 
tightening of the gradient can be seen over the western val- 
ley wall when compared to the eastern wall, but only over 
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tha northern half of the valley. At tha highest two levels, 
gradients over the two valley walls are nearly identical so 
that the apparent changes in tamperatura are probably due to 
the effect of terrain variation on tha slope of the sigma 
surface. 

The cold air continues to penetrate southward for the 
next twelve hours reaching its maximum southward extent at 
108 h (Fig. 17). At 24°N the minimum tamperature is 298K, a 
-12K change frcm 84 h. The cold air does not penetrate much 
beyond the southern valley border. At 20°N the temperature 
change is approximately -2K in 24 h. i^hile southward move- 
ment of the cold air ceases at 108 h, it doss expand east- 
ward, apparently in response to the eastward movement of the 
anticyclone located north of the vallay. The cold air 
tongue does not dissipate as rapidly as the pressure ridge, 
and a pool of relatively cool air remains in the western 
valley region until the arrival of the next mid-latitude 
cyclone. 

At the higher levels, sigma = 7/12 and 9/12, a small 
tongue of cold air forms between 78 and 34 h above the west- 
ern valley wall. The southward extant is slightly less than 
that at the lowest level. As the low-level cold air tongue 
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moves rapidly southward aftar 84 h, the cold air at rhese 
levels moves very slowly, only four degrees between 94 and 
102 h at sigma = 7/12. Little, if any, eastward expansion 
is seen. Beyond 108 h the very slight gradients are diffi- 
cult to discern, primarily due to the increasing small scale 
noise. At sigma = 9/12 the cold air more closely follows 
the behavior of the cold air tongue at the lowest level 
showing continued south and eastward expansion on the same 
time scale, although of less intensity. 

To verify that the changes in the potential temperature 
fields at the upper-levels are essentially due to the use of 
the sigma coordinate system, the temperature (T) fields were 
analyzed at 850, 700 and 500 mb between 84 and 102 h. At 
850 mb the intrusion of ccld air can easily be traced. At 
700 mb the cold air tongue is not discernable although this 
may partially be due to the adiabatic warming effects appar- 
ent at both the 700 and 850 mb levels. At 500 mb some small 
scale noise is introduced by flow over and around the valley 
but there is no evidence of the southward movement of cold 
air. Therefore it is reasonable to assume that the changes 
in the potential temperature fields at the higher levels are 
due to the effect terrain variation has on the slope of the 
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sigma surface. The apparent movement of colder potential 
temperatures southward, then, is merely the result of the 
sigma surface being at a lower elevation in the valley 
region. To seme extent this effect is seen near the eastern 
valley wall at the lower levels as well. 

At 78 h, when the valley is only partially completed, a 
southward deflection of the northwesterly winds ahead of the 
ridge is readily apparent. By 84 h the valley is complete 
and the north winds are rapidly moving southward (Fig. 18). 
Along the valley wall the leading edge of the wind surge is 
highly ageostrophic . A strong southwesterly flow has set up 
in the eastern region of the valley presumably in response 
to the tight pressure gradient that exists there. It is 
interesting tc note the small cyclone that has developed 
downwind of the western valley wall. It appears to be 
slightly better defined than the pre-existing low to the 
northeast and may result from the combined effects of the 
shear line and of cyclonic flow induced in the descending 
air by the conservation of potential vorticity. The cyclone 
dissipates as the north wind continues to push south and 
eastward into the valley. The shear line remains well 
defined between the opposing northerly and southwesterly 
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flows. The strong, cross-isobaric southwest winds in the 
eastern valley coincide with a very tight pressure gradient 
and may represent a sloping front in the valley. 

At 96 h the north wind regime extends rhe length or the 
western valley wall and is spreading eastward across the 
valley (Fig. 19). At 102 h t he north winds have reached 
their maximum extent; the wind regime occupies fully half 
the valley but does not extend south of it. The flow along 
the valley wall is no longer age ostrophic. The north winds 
entering the valley between 84 and 102 a are identifiable 
only at the very lowest levels, in agreement with the defi- 
nition of a monsoon surge as a low-level phenomenon. By 108 
h the anticyclone has moved sufficiently eastward that the 
light northerly winds associated with the ridge line have 
replaced the strong northwest flow previously located north 
of the valley (Fig. 20). The flow into the valley is 
greatly reduced and, without this support, the surge weakens 
and dissipates. The eastern region continues to be domi- 
nated by southwest flow. 

It is clear that the movement of the north winds through 
the valley correlates well with the development of the cold 
air tongue and pressure ridge discussed earlier. 
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In suaunary, savsral observations can be made from rhe 
analysis of these fields. First, all three fields shew the 
very rapid southward movement of a cold, high pressure air-- 
mass along the western valley wall between 84 and 96 h. 
Much slower southward movement continues through 108 h. 
Second, the airmass expands eastward across -he valley in 
response to the continued eastward movemeni: of zhe anticy- 
clone. And third, the airmass never penetrates the tropical 
regions but remains within the confines of the valley. 

C. EXPERIMENT 4 

In this experiment the mountain range is again centered 
at 30°N. The valley, built between 72 and 84 h, is 24 
grid-lengths wide, ie. two-thirds of the mountain range in 
length. The terrain then is actually the small east-west 
oriented mountain range displayed in Fig. 21. As in Experi- 
ment #3, the critical period for the surge is between 82 and 
108 h while the anticyclone is located northwest of the 
valley. 

The first 72 h of development are described by Control 
Run #2. As in the previous experiment the creation of the 
valley causes an oscillation in the pressure (ie. mass) 
field. The oscillation in this experiment is larger than 
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rhat of the previous case since the reduction of terrain is 
greater. The rapid increase in cyclogsnesis also occurs 
earlier, supporting the assumption that the mountains exert 
a stabilizing influence on the disturbanos. 

At 84 h a small southward extension of the high pressure 
ridge can be seen over the eastern mountain slope (western 
valley wall). A low is centered over the northern valley 
with a weak trough extending southwest over the valley wall. 
A tight gradient develops south of the low, extending over 
the eastern valley wall and continuing along the northern 
mountain slopes. The trough is strongest at 78 h and is 
subsequently pushed southward by the expanding high pressure 
ridge. The ridge advances southward during the next twelve 
hours reaching the southern valley edge at 96 h (Fig. 22). 
Thus far movement of the ridge is quite similar to Experi- 
ment #3, although it does appear to be slightly weaker and 
not as readily discernable. 

Between 96 and 108 h the mid-latitude ridge weakens con- 
siderably. Concurrently, a strong pressure gradient builds 
along the western edge of the ridge, in part due to a trough 
extending southward from the intense cyclone located to the 
northwest (Fig. 23). A similar gradient is seen during the 
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same period in Experiment #3 but is nor nearly as intense. 



At both 102 and 108 h there is evidence that th^ 


: ridge 


extends slightly beyond the edge of the valley. 





The development and movement of the ridge agrees very 
closely with that of the earlier experimenr using a smaller 
valley. In both cases the ridge begins development while 
the valley is growing, moves southward along the western 
wall at the same rate, reaches the southern edge at the 96 h 
point and weakens when the anticyclone moves to a location 
due north of the valley. There is, however, a slight dif- 



ference in intensity. The third experiment shows 


a well 


defined, easily discernable ridge from 34 to 108 h. 


In the 


present experiment the gradient along the western 


side of 


the ridge is greater, but the ridge itself is not 


as well 


defined. 





A weak tongue of cold air can be seen along the western 
valley wall at 84 h (Fig. 24). An increase in the tempera- 
ture gradient is apparent in the extreme western valley at 
the lower three levels although it is very weak at the 
uppermost of these (sigma = 7/12). These levels correspond 
to pressure heights lower than 600 mbs in a standard atmos- 
phere. Through the next twelve hours the cold air tongue 
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moves southward and strengthens. 3y 96 h the cold air has 
reached the southern edge of the valley (Fig. 25). The cold 
air tongue is much better defined than the pressure ridge. 
It is nearly identical to the cold air tongue seen in the 
western valley for Experiment #3. The minimum temperature 

at 24®N has fallen from 310K at 84 h ro 303K at 96 h. The 
temperature change is primarily limited to -he western third 
of the valley, temperatures in the eastern cwo-thirds have 

4 

remained relatively constant following completion of the 
valley. 

The cold air continues to move south and eastward 
through the next twelve hours. Maximum southward extent is 
reached at 108 h. 4t 24®N the minimum temperature is 298K 
and represents the same change in temperature, -12K, as was 
seen in Experiment #3 (Fig. 26) . 'rfhile southward movement 
of the cold air has ceased, it does continue to expand east- 
ward well past this time. The primary axis of the cold air 
tongue remains along the western valley wall and retains its 
identity through 120 h. A pool of cool air remains in the 
western valley for another twelve hours. At sigma = 9/12 
the cold air is discernabla well past 108 h. At higher lev- 
els the small scale noise has more effect and the cold air 
is difficult to trace past this time. 
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The upper-levels display the same behavior seen in 
Experiment #3. A slight increase in eastward expansion at 
these levels is seen, possibly due zo the much wider valley. 

As with the pressure and temperature fields the develop- 
ment of the north wind surge in this experiment closely fol- 
lows that of Experiment #3. At 84 h the norrh winds have 
reached mid-valley and continue to spread rapidly southward. 
Again the leading edge of the wind surge is ageost rophic. A 
strong southwesterly flow has set up in the valley appar- 
ently due to the tight pressure gradient located there (Fig. 
27) . The small cyclones seen in Experiment #3 are not 
apparent here although cyclonic flow exists, possibly due to 
the convergence of the southwest flow from the tropics and 
the opposing flow out of the mid-latitude anticylones. This 
cyclonic flow weakens considerably after 34 h. 

At 96 h the north wind regime extends almost the entire 
length of the valley and is spreading eastward. The south- 
ward extent of the north winds along the valley wall is 
identical to that of Sxperiemtn #3 (Fig. 28). Further east 
the perimeter of the north wind regime is nearly two degrees 
farther south of that reached in the earlier experiment. 
The mid-latitude ridge weakens during the next six hours. 
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and with its ccntinued eastward movement, results in a 
decrease in the northwesterly winds supporting the surge. 
Thus only a very little southward movement takas place 
between 96 and 102 h, the time of maximum southward extent 
of the north winds. ' The ageostrophic component of the flow 
weakens rapidly as southward movement along the valley wall 
slows. The winds become much weaker but do not decrease in 
coverage through 108 h (Fig. 29). k slight eastward move- 
ment is seen as the ridge moves past the mountains and north 
of the valley. After 108 h the north winds diminish and 
move eastward in response to the weakening and continued 
eastward movement of the ridge. After 120 h the wind surge 
is no longer discernable. 

The surge is identifiable only at the two lowest levels. 
At sigma = 9/12 the surge is weaker and lass extensive than 
at the lowest level. It dissipates very rapidly after 102 h 
as the ridge at this level moves north of the valley. 

In this experiment a slight increase in the eastward 
extent of the north winds is noted when oompared to Experi- 
ment #3. Since topography is the only parameter changed, it 
is reasonable to assume the increase in eastward extent is a 
result of the expanded valley. It should be noted though 



45 



that while the valley size has double!, ohe corresponding 
increase in longitudinal coverage by the north winds is 
extremely small. 

Several observations can be made from the results of 
Experiment #4. First, all three fields show the rapid 
southward movement of the cold, high pressure airmass along 
the western valley wall between 84 and 96 h. Subsequent 
southward expansion is much slower. Second, the airmass 
expands eastward across the valley in response to the east- 
ward movement of the anticyclone. And third, there is a 
tendancy towards cooling and rising pressures just south of 
the valley although the north winds remain within the valley 
confines. 

D. EXPERIMENT 5 

In this experiment the mountain range is centered at 
38®N, eight degrees farther north than in previous experi- 
ments. The valley is twelve gridlengths wide, as in Experi- 
ment #3. Terrain is shewn in Fig. 33. Cyclogenesis is 
considerably slower than that seen in either of Experiments 
3 or 4. This is due to the placement of the northern moun- 
tain slopes at the same latitude as the maximum amplitude of 
the disturbance. Growth of the perturbation is prevented by 
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effscrively blockiag the north-south heat flux. 



For this 



reason the initiation of valley growth was delayed six hours 
to allow additional time for development of the disturbance. 
The growth period for the valley is now between 78 and 90 h. 
As with the previous .experiments the creation of the valley 
causes an oscillation in the pressure field, the approximate 
amplitude of that in Experiment #3. 

At 90 h a strong low is located in che northeast of the 
valley trough extending southwest over the valley wall. 
This results in a tight pressure gradient in the eastern 
valley and along the northern mountain slopes similar to the 
gradients found in previous experiments. The mid-latitude 
anticyclone is neither as strong nor as well developed as in 
Experiments 3 and h, and although there does appear to be a 
slight southward extension of the high pressure along the 
western valley wall, it is very weak. During the next 
twelve hours a ridge of high pressure builds along the val- 
ley wall, however it is difficult to determine whether the 
high pressure originates in the anticyclone to the northwest 
or in the tropical regions south of the valley (Fig. 31). 
As the analysis of the wind fields will later show, the high 
pressure air does originate in the mid-latitude anticyclone. 
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At 102 h the anticyclone is 



ili-definei and rha aounzain 



region to its south populated by small high and low centers. 

Between 102 and 1 14 h the ridge moves east into the val- 
ley and strengthens, this is par ticularily noticeable at 114 
h (Fig. 32) . It appears the anticyclone has slipped 
beneath the cyclone and is re-forming in the northarn val- 
ley. During this time a weak zrough can be seen forming 
over the northwestern valley wall. 

At 120 h the anticyclone is locazad in the northeast 
valley region wizh a weak ridge extending southwest througn 
the valley and a stronger, better defined ridge to the 
northwest. Low pressures are found over the weszern valley 
wall due to a weak trough forming ahead of the mid-latitude 
cyclone. 

In earlier experiments the high pressure ridge developed 
and strengthened along the western valley wall only to dis- 
sipate as the anticyclone moved to a position northeast of 
the ridge. In this experiment the ridge is weaker, not 
well-defined and moves eastward with the anticyclone. Since 
the northern valley is located at the latitude of maximum 
amplitude for the disturbance, creating the valley allows 
the anticyclone to reform in this area as the weak ridge 



moves eastward into the valley. The same reasoning can be 
used to explain the strong development experienced by the 
low during the valley growth period; it was locaned along 
the northern mountain slopes at the same longitude as the 
valley. 

A weak tongue of cold air can be seen along the western 
valley wall at 90 h (Fig. 33) . This gradient is found in 
the lowest three levels, weakening with height. At the 
upper-levels (sigma = 5/12, 3/12, 1/12) nearly identical 

gradients exist at both ends of the valley and are due to 
the effects of terrain on the slope of sigma surfaces. The 
cold air tongue strengthens and moves rapidly southward in 
the following six hours. By 102 h the cold air has reached 
the southern edge of the valley (Fig. 3’4» . The development 
and movement of the cold air tongue is much better defined 
than that of the pressure ridge. It is nearly identical in 
development and intensity to the cold air tongues seen in 
Sxperiments 3 and 4. The minimum temperature at 320N (a 
comparable latitude to 24°N relative to the position of the 
mountains) is 295K compared to 306K at 90 h. This is a 
temperature change of -11K compared to -12K and -10K for 
Experiments 3 and 4 respectively. Again the temperature 
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change is primarily limited to the wsstarn third of the 
valley. 

The cold air continues to move south and eastward during 
the next twelve hours, reaching its maximum southward extent 
at 1 14 h (Fig. 35) ► At 32®N the minimum zemperarurs is 
293K, a temperature change of -13K compared to -12K in 
Experiments 3 and 4 for a 24-h period following complezion 
of the valley. Eastward expansion of the cold air continues 
past 114 h. In the previous experiments the axis of the 
cold air tongue remained along the western valley wall. in 
this experiment the cold air moves eastward as well, parti- 
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The upper-levels display the same behavior seen in 
Experiments 3 and 4. A very slight movement of the northern 



portion 


of the cold air axis 


is seen 


at sigma = 9/12, but 


no 
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movement 


is seen at 
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t»mp«5ratur € fields difficult past 114 h. The highest levels 
(sigma = 5/12 to 1/12) show ao sign of southward movement of 
cold air at any time. 

As in the previous experiments southward deflection of 
the northwest winds into the valley occurs before the valley 
is completed. By 90 h the north winds have reached mid-val- 
ley with a shear line developing over the valley floor in 
advance of the northerly flow. North winds along the valley 
wall are strongly ageost rophic. Again a strong southwest 
flow has set up in response to the tighn pressure gradient 
in the eastern valley. Although the valley is the same as 
that used in Experiment #3, the subsequent pattern of devel- 
opment is slightly different. At 96 h the northwest winds 
reach the eastern valley wall, whereas, in Experiment the 
strong southwest flow in this region prevents the northwest 
winds from spanning the valley. The wind surge is actually 
strongest at this time, and has its greatest southward 
extent (Fig. 36) . By 102 h the ridge is north of the valley 
wall and the winds are beginning to weaken, the flow is no 
longer ageostrop hie. The nornh winds never reach the south- 
ern edge of the valley (Fig. 37). The surge continues to 
weaken rapidly and at 108 h the anticyclone is reforming 
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wirhin the valley. North winds exist in the central valley 
region but dc not coincide with the positon of rhe cold, 
high pressure ridge. Six hours later rhe anticyclone is 
esrablished over the north central valley and no evidence of 
the surge exists (Fig. 38). It is clear that positioning 
the mountain range at 38°N has a detrimental effect on the 
development of the surge. 

The surge is identifiable at only the two lowest levels. 
At sigma = 9/12 the surge is less extensive and weaker than 
at the lowest level. At sigma = 7/12 and 5/12 a weak 

cyclonic turning of the winds is apparent but no definite 
southward flow can be seen. The cyclonic flow disappears as 
the anticyclone reforms over the valley. 

In summary the following observations can be made from 
this experiment. First, all three fields show the rapid 
southward movement of cold, high pressure air along the 
western valley wall during the first six hours following 
valley completion. Second, the movement and development of 
the cold air tongue is much better defined than that of the 
pressure ridge and wind surge, both of which weaken consid- 
erably after 102 h. Third, the pressure ridge and wind 
surge axis move eastward off the valley wall due to the 
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development of the anticyclone within the valley, only the 
northern most portion of the cold air tongue moves eastward. 
And fourth, although there is a slight cooling and pressure 
rise south of the valley, the north winds remain well within 
the valley. 

E. EXPERIMENTS 6 AND 7 

The valley was included as part of the initial terrain 
configuration for these two experiments. The mountains are 
centered at 30<^N and the valley is eight gridlengths wide, 
approximately one-guarter the length of the mountain range. 
This is the original mountain -valley configuration used with 
the model, the larger valley was used in subsequent runs to 
avoid interference between the strong southwesterly flow in 
the eastern valley and developing southward flow in the 
west. The southward movement of a cold, high pressure air- 
mass occurs with both valley sizes, however the airmass is 
more easily seen, and so more easily analyzed, with the 
larger valley. 

Only the mean flow is initialized in Experiment #6. The 
purpose of this run is to confirm the assumption that the 
monsoon surge, if it occurs, is due to baroclinic wave and 
not a result of the interaction of the mean flow with the 
terrain. 
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II 



As expected no "surge" occurs. The current rhar does 
form in the valley is southwesterly, flowing from high to 
low pressures. The strong southwesterly flow seen in Exper- 
iments 3 through 5 is in part due to the mean pressure gra- 
dient that exists between the northern and southern valley 
borders. By 30 h a small cyclone has formed over the wesr- 
ern valley wall, probably due to rhe conservation of poten- 
tial vorticity. This indicates that the small cyclone seen 
here and in ether experiments is a result of the mean flow 
interacting with the terrain. 

Experiment #7 is initialized with the perturbation. 
During the first portion of the experiment, while the dis- 
turbance is weak, southwesterly flow forms in the valley in 
response to the existing north-south pressure gradient. 
Again at 30 h a small cyclone has formed over the western 
valley wall. During the next twelve hours the center of the 
cyclone moves ever the valley floor althoagn the north winds 
remain located over the valley wall. The mid-latitude dis- 
turbance continues to strengthen and by 56 h the northerly 
winds are being enhanced by the northwest flow ahead of the 
approaching ridge. The northerly flow strengthens as the 
ridge moves eastward, at 96 h it extends slightly beyond the 
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valley edge, however this is its maximnn southward extent. 
Subsequently the mid- la ti tude ridge axis moves north of -he 
valley and the supporting northwest flow diminishes. The 
flow along the valley wall is never ageostrophic during this 
time. 

The cyclone is readily apparent in the sea-level pres- 
sure fields. However, not until the ridge approaches the 
valley does a high pressure ridge build along the valley 
wall (72 h) . The ridge weakens as the anticyclone moves -o 
a position north of the valley. 

The potential temperature field shows very little cool- 
ing along the western valley wall compared to previous 
experiments. Minimum temperature along 2'4®N an 90 h is 299K 
compared to 303K at 66 h, a decrease of only 4K in 24 h. A 
very small temperature change when compared to Experiments 3 
through 5. 
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Figure 6. Potential TeniDeratura , 
Sigma = 11/12Z 
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59 







Figurs 7. Model Terrain, Expert aient 2. 
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Figure 9. Sea-Levsl Pressure, 84 and 95 h, Sxp. 2. 
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Figure 10, 



Potential Tesiperature^ 
Sigma = 11/12. 
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Potential Temperature, 
Sigma = 11/12. 



34 and 95 h, 2xp. 2. 
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Figure 16, Porential Temperature, 96 h, Sxp. 3, 
9/12, 11/12. 
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Figure 17. Potential Temperature, 108 h, Exp. 3. 
9/12, 11/12. 
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Figure 21. Mcdel Terrain, Experiiuenr 4. 
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Figure 23, Sea-Level Pressure, 108 h, 3xp* 4* 
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Fiqure 24. Potential Temperature/ 34 h, 3xo. 4. 
9/12, 11/12. 



Sigma = 



77 



90 -I 



67 - 



414 H 



21 “ 



315 



-2 - 




rso 



-67 






-21 



r --2 



23 



45 



66 



I 

90 



90 -T 



67 - 



44 - 



21 - 



-2 - 




265.0 

270.0 



275.0 



280.0 

285.0 

230.0 




310.0 



305.0 



r 90 



r 57 



-44 



-21 



--2 



Figure 25. 
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9/12, 11/12. 
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Figure 28, 
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Figure 29 . 



Wind Vectors, 108 h, SxDerioienr 4, Siqma 
9/12, 11/12. 
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Figure 30. Model Terrain, Experiment 5. 
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Figure 33. Potential Temperature, 90 H, Exp. 5. Sigma = 
9/12, 11/12. 
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Figure 34. Potential Teaperarure , 102 li, 3xp. 5. 
9 / 12 , 11 / 12 . 
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7. 



CONCLO SIGNS 



The interact ioa of a mid-latitude baroclinic wave wirh 
several different mountain-valley configurations vas exam- 
ined. The purpose was to determine the contribution made by 
a well-developed baroclinic wave to the initiation of mon- 
soon surges. It was found that these waves can initiate a 
monscon surge although such surges are wealc and limited in 
extent. 

The mountains in each experiment are an east-west range 
located south of the mid-latitude wave. As was first noted 

in Experiment #2 the mountains appear to have a stabilizing 
effect on the growth of the wave. The effect is especially 
evident when the mountains are placed farther north, nearer 
the center of the disturbance. The reduced growth rate is 
probably caused by the blocking of the v-component of the 
wind by the mountain barrier to the south. This would slow 
the wave growth because the energy conversion is porpor- 
tional to V'9*. Walker (1982) also found this effect with a 
disturbance moving along the north side of a mountain range. 

Previously, several observations were made about the 
nature of monsoon surges. These include: (1) the occurrence 

of rapid, ageostrophic flow from high to low pressure, (2) 
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strengthening of the north-northeasterly winds, (3) a sharp 
drop in surface teaperature and (4) rhe shallow structure of 
the surge. To seme extent, each of these cnaracterisr ics is 
simulated by one or more of the experiments. 

Experiment #6, in which only the zonal mean flow inter- 
acts with the terrain, shews no evidence of a surge. All 
other experiments, in which the perturbation is added to the 
mean flow, do produce weak monsoon surges. This indicates 
the eastward- moving anticyclone associated with the mid-la- 
titude baroclinic wave is the source of the surges occurring 
in this study. 

In Experiments 3 and 4 the mountains are centered at 
30°N, the western edge of each valley is located at the same 
point although the fourth valley is twice the size of the 
third. In each case a cold, high pressure airmass moves 
southward along the western valley wall during the twelve 
hours following completion of the valley. The leading edge 
of the surge is ageostrophic. 

Figs. 39 and 40 show the change of potential temperature 
with time in the valley for Experiments 3 and 4. The pro- 
files extend from 140N to 46°N along the western valley 
floor (55° long) , the valley is located between 22°N and 
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38°N. In each case little, if any, change in c emper ac uc e 
occurs at 14°N although at 18°N the temperature decreases 7 k 
between 96 and 108 h. In the valley the maximum temperature 
change occurs earlier, between 84 and 96 h, and is much 
larger, -16K and -17K at 26®N. The surge is well defined 
within the valley but dissipates rapidly south of 22°N. 
This is indicative of the weak character of the monsoon 
surges in this study since the cold air in observed cases 
has been traced well into the tropical regions while in 
these experiments it cannot be found even eight degrees 
south of the valley. 

As can be seen the results of the two experiments are 
nearly identical, indicating the size of the valley is not 
an important consideration. 

The oscillation in the pressure (mass) field is appar- 
ently a gravity wave, the result of an inward rush of air 
filling the "vacuum” produced by the creation of the valley. 
It has roughly a 24-h period, which may be related to the 
twelve hour period in which the valley is built, and it 
slowly dampens during the subsequent 60 h of the experiment. 
It has no observable effect on the pressure gradients and 
was not evident in either the potential temperature or the 
wind fields. The oscillation also occurs in Experiment #5. 
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In Experiment #5 the rapid southward movement of the 
cold airmass occurs only during the first six hours after 
valley completion. This is due zo the placemenr of rhe 
mountains at 38°N resulting in a much weaker disturbance. 
Fig. 41 shows the change in potential remperarure with time 
in the valley. The valley is located between 30°N and 46°N, 
the profile is along the western valley floor. 

The decrease in temperature is evident in the valley 
between 90 and 102 h. As in Experiments 3 and 4 the cold 

air can be traced only four degrees south of the valley. 
Unlike the two previous experiments warming occurs during 
the second twelve hour period, due to the re-development of 
the anticyclone within the valley. The southerly winds in 
the western region of the anticyclone prevent the flow of 
cold air into the valley and effectively "shut off" the 
surge. 

The monsoon surges initiated in these experiments are 
evident only in the lowest sigma levels, at sigma = 11/12, 
9/12 and 7/12. In a standard atmosphere, over flat terrain, 
these correspcnd to pressure heights less than 675 mb. The 
surge is strongest at the lowest level and dissipates with 
height, becoming very weak at sigma = 7/12. This is 
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observed in both the wind ana potential temperature fields 
in all experiments and agrees closely wirh the definition of 
the monsoon surge as a low-level, or surface, phenomenon. 

These results suggest that the mid-latitude baroclinic 
wave can initiate monsoon surges. The surges are weak and 
limited to the valley but otherwise exhibit the characteris- 
tics of observed monsoon surges. The weakness of these 
surges indicates other forcing mechanisms are necessary in 
order to adequately model the monsoon surge. These may 
include a strong, stationary high located northwest of the 
valley similar to the anticyclone associated with the winter 
monscon. Other analytical studies indicate that an upper- 
tropospheric, planetary scale wave number three is necessary 
for the occurrence of strong monsoon surges. Such a wave 
positioned with the ridge to the west of the surface anticy- 
clone would result in strong northwesterly flow over the 
high pressure center. This would initiate cold air advec- 
tion throughout the troposphere and increase subsidence over 
the surface anticyclone leading to stronger surges. Weak 
monsoon surges were associated with near zonal flow at 
higher altitudes, this describes the upper-level flow used 
in these experiments. The inclusion of one or both of these 
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plienomena may result in stronger, more extensive 
surges than these observed in this study. 



monsoon 
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gure 39. Potential Tempera-cure vs Larituie, Exp. 

Vertical lines reoresent valley boundari 
84 h, (b) 96 h and (c) 108 b. 
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Figure UO. Potential remoerature vs Latitude, 2xd . 

Vertical lines represent valley boundari 
84 h, (b) 96 h and (c) 108 h. 
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Figure 41. Poterrial Temperature vs Latitude, Exp. 5. 

Vertical lines represent valley boundaries. (a) 
90 h, (b) 102 h and (c) 14 a. 
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APPENDIX A 



aaOEL PRIMITIVE EQUATIONS 

Equation of Stats: 
pa = RT 



Hydrostatic Equation: 



ii _ 

da ~ 



-ira 



Horizontal Momentum Equations: 

n/dv\ + fk nv + 7 (n$) - vn §_(^a) = o 

\dt/ ^ Sa 

Continuity Equation: 

5n + 7 • (nv) + ^(na) = o 

5t ^ So 
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1st Law of Thermodynamics: 



s (ncpine) + V* (nvcpine) + §_ 

6 5a 



(IlaCplnS) 
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LIST OF SYMBOLS 

radius of the earth 
amplitude of the perturbation 
specific heat at constant pressure 
ccriclis parameter 

value of Coriolis parameter ax 45N 
force 

gravitarional acceleration 
wave number 

pressure (P = P + P' ) 

surface pressure at U5N, 1013.25 mb 

maximum model pressure height, 125 mb 

gas constant for air 

temperature 

surface xamperature at 45N, 288K 

zonal velocity component (U = U + O') 

maximum velocity, surface: 5 m/s, upper: 65 m/s 

meridianal velocity component (V = V + V') 

terrain height 

specific volume 

halfwidth of jet 

dry adiabatic lapse rate, 6.5 X/!cm 

potential temperature 

longixude 

3.14159 

normalized pressure 
vertical coordinate 

vertical velocity in sigma ccoridinate system 
model time 

model time at which terrain growth begins 
model time at which terrain reaches final heigh 
latitude 
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<i>j latitude of jet 

latitude at which mountain range is 
geopotential, also current terrain 
terrain height, 3000 m 
velocity of the earth 






maximum 

angular 
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